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The results of G. Su and M. Suzuki [Phys. Rev. B 54, 8291 (1996); ibidem 57, 13367 (1998)] for 
the spin and pseudo-spin symmetry properties of the Hubbard model are reexamined. We point out 
that the exact relations they have found are valid down to zero temperature and that their solutions 
for both spin and pseudo-spin correlation functions are incorrect. 
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In a recent paper G. Su and M. SuzukiB analyze the 
spin symmetry properties enjoyed by the Hubbard model 
in presence of an homogeneous external magnetic field. 
They derive at finite temperatures an exact relation con- 
necting the spin correlation function to the magnetiza- 
tion. Also, the authors claim to have found without any 
a priori assumption at least one of the exact solutions for 
the magnetization as function^ of the applied held. This 
result follows previous worksaQ for the pseudo-spin coun- 
terpart where a solution for the pseudo-spin correlation 
as a function of the filling has been assessed. 

In this Comment we clarify the issue of applicability 
of the exact relations, by use of the equation of motion, 
showing that they are valid also at zero temperature. 
Furthermore, we show that the pretended solutions for 
both spin and pseudo-spin correlation functions are in- 
correct. 

The Hubbard model in presence of an homogeneous 
external magnetic field h reads as 

H = J2 (*« - A*%)c t (i)c(j) + Uj2 n l(i)ni(i) 

ij i 

(1) 

i 

where s z (i) is the third component of the spin density 
operator. 

Let us introduce the total spin operators 
S+ = X)4(i)q(i) 

i 

s~ = J2°\(i)°r® (2) 

i 

s z = -^2 ( c ! W c t w - 4 W c i (*)) 

i 

and the thermal retarded Green's function 
S+- (t-f) ={K[S+(t)S- (*')]) 

= — [ lioe- l ^ t - t ')s + -{uj) (3) 

2^ J-oo 

By means of the Hamiltonian (|l|) the spin operators sat- 
isfy the Heisenberg equations 



i — S ± = ±2hS ± i — S z = 
dt dt 



Then, we have 



— S + ~(uj) = 2m 

N [UJ> u;-2h + it] 



(4) 



(5) 



where N is the number of sites and m is the magnetiza- 
tion per site 



(6) 



In presence of an external magnetic field the spin sym- 
metry is explicitly broken, [cfr. Eq. (Q)], and the-propaga- 
tor (lj) exhibits a massive collective modeel u = 2h. 
When h = and m ^ the collective mode becomes 
gapless, in accordance with the Goldstone theorem. 

From Eq. (||), by standard methods, we obtain the spin 
correlation function 

(S + (t)S-(t'))= 2m i e _ p _ 2 , h (7) 



N 



where (3 = 1/k^T. Similarly, we derive 

1 2me 2ih ( t - t ') 
jj{S-{t)S+{t>)) = - 2m i e _ eWh (8) 

Let us note that (M) and (||) satisfy the KMS relation: 
(S+(t)S-(t' + i f3)) = (S-(t')S+(t)) (9) 
In the static case Eqs. (^) and (@) become 



— (S + S-) = m [coth (/3h) + 1] 

— (S~S + ) = m [coth (flh) - 1] 



(10) 



These are the exact relations derived by Su and SuzukiB; 
however they are not restricted at finite temperature, but 
hold also at T = 0. In particular, for finite magnetic field 



lim —(S+S-) = 2r 
lim — (SS+) = 



(11) 
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FIG. 1. Magnetization m as a function of filling n, applied 
magnetic field h and temperature T. 

On the basis of the relations (|l(]) Su and Suzuki0 pro- 
mote as one of the exact solutions for the magnetization 
as a function of the applied magnetic field the following 
expression 



m (h, T) = - tanh (fih) 



(12) 



where n is the particle density. Also, they stress that 
other solutions, if they exist, might have similar forms. 
At the opposite, the solution depicted in (1^) is clearly 
wrong except for the limiting case of half-filling and infi- 
nite U. Indeed, the pretended solution (|l^) can be falsi- 
fied by looking at two exactly solvable limits of the Hub- 
bard model. That is, the noninteracting [i.e. U = 0] and 
atomic [i.e. t^j — 0] ones. 



Non interacting case 



It is direct to see that 



n = 1 -^h I ddfc PXk)+n(k)] 

2 (2ir) Jn B 



m — — tanh (f3h) 



1 - 



it 



(2^) a Jn 



d d feT T (k)T;(k) 



(13) 
(14) 



51 is the volume of the unit cell, d is the dimensionality 
of the system and Ob is the first Brillouin zone. We put 
T CT (k) = tanh [[3 E a { k)/2] with energy spectra 



£ T (k) = -fi - 4ta(k) - h 
EiQs) = -fx - 4to(k) + h 



(15) 



where a(k) = l/dz~2i = i cos (^« a )' ^ i s ^ ne hopping integral 
and a is the lattice constant. 

In the two-dimensional case the magnetization is 
shown (cfr. Fig. |l|) for different values of the parame- 
ters n, T and h. The solution (p"2|), proposed in Ref. [I], 
obviously is not a solution for the case U = 0. 



Atomic limit 

In this case it is easy to show that 

_ 1 - \ (Ti +T 3 ) + \ (1- T 3 ) (Ti - T a ) 



2m 



1 - | (Ti - T 2 ) (T 3 - T 4 ) 
I (1 - Tj) (T 3 - T 4 ) 
1 - i (T x - T 2 ) (T 3 - T 4 ) 
| (T 3 - TQ + 1 (1 - r 3 ) (Ti - r 2 ) 
1 - \ (T x - T 2 ) (T 3 - T 4 ) 
_ Kl-TO^-^) 
1 - i (T x - T 2 ) (T 3 - T 4 ) 

where 1$ = tanh ((3 Ei/2) with energy spectra 

i?i = — /i — /i 
£ 2 = -H + U - h 
E 3 = -fi + h 
E 4 = -n + U + h 



(16) 



(17) 



(18) 



We note that for /i = U/2 the previous expressions be- 
come 



n = 1 



Ti +T 2 



= Thjl + Ty) 

4 [! + 1 (Ti-T 2 )] 2{l+TuTl) 



with 



T v = tanh (P U/4) 7), = tanh (/3 h/2) 



(19) 



(20) 



In particular, in the limit of large U and for finite tem- 
perature 



l + T, 



, — — tanh (/? h) 



(21) 



in agreement with the result of Ref. [|. 

An intrinsic symmetry of the Hubbard model is the 
pseudo-spin SU{2) symmetrjfl that combined with the 
spin SU{2) one yields the [SU{2) <g> SU{2)]/Z 2 = 5*0(4) 
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symmetry group. The generators of this transformation 
are given by the total pseudo-spin operators 



P+ = £ e ;Q-R 4c ; (i)c t w 

i 

p-=£>-^ Ci (i)c T (i) 

i 



(22) 



where Q = (it, it). These operators satisfy the Heisenberg 
equations 



d 



i—P ± = ±(2 M -E/)P ± 



*^P 

l dt z 



(23) 

Let us consider the thermal retarded Green's function 
P+-(t-t') =(K[P + (t)p- (*')]) 

/+oo 
dwe"^( t - t ')p+-( w ) (24) 



(25) 



By means of the equation of motion (|23|) we find 

1 _| n — 1 

iV ^ = w - (2/x - C/) + 



The Hamiltonian has pseudo-spin SU(2) symmetry only 
at half-filling; when rt ^ 1 the symmetry is explicitly 
broken. Again, this is seen in the propagator .where a 
massive collective mode u> = 2/i — U is observed! From 
(Bq) we obtain the correlation function 



-(P+(t)p- (t'))= { 



and similarly 
1 



N 



(p-(t)p+(tO) = - 1 1 ) _ e ^. u) 



(26) 



(27) 



It is easy to see that (|26|) and (|27j) satisfy the -075 re- 
lation. In the static case: 

1 (P+P-) = \{n - 1) [coth {0(2n - C/)/2) + 1] 

1 (P-P+) = i(n - 1) [coth (/3(2 M - EO/2) - 1] (28) 



which give 



(P+p-) = e^- c/ ) (P-P+) 



(29) 



These exact results relates the pseudo-spin correlation 
functions to the particle number n and are a manifesta- 
tion of the intrinsic symmetry. These relations generalize 
at T = the results previously obtained by Sutj. 

Under the particle-hole transformation we have the fol- 
lowing relations 



u 




FIG. 2. C as a function of filling n, intrasite Coulomb re- 
pulsion U and temperature T. 



fj,(2-n) = U-n(n) C + ~(2 - n) = C _+ (n) (30) 

where we put C+~ = (P+P"), C~+ = (P~P + ). Use of 
Eq. @ leads to 



C+"(2 -n) = e-^W-^C+'Cn) 



(31) 



By making use of the transformation properties ( j30| ) , it 
is easyto see that the expression (^8|) satisfies the prop- 
erty (|3l]). Actually, this is a manifestation of the inti- 
mate interrelation between pseudo-spin and particle-hole 
symmetries! 

In Ref. H the particle-hole symmetry is tautologically 
used as a supplementary equation and the following solu- 
tion for the pseudo-spin correlation function is presented 
for the case h = 



N 



(P+P-) 



C(T) 



(32) 



where C(T) is an unknown function of temperature only. 
When ( p2] ) is used in ( p8f ) one is lead to the following 
equation for the chemical potential 



n = 1 + C(T) tanh [/3(/i - 172)] 



(33) 



This equation is incorrect. For example, let us consider 
the limit of small temperature. Then, ( |33| ) would give 
n — > 1 — C(0), which is clearly wrong. 
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In Fig. H we present the function C(n, T, U) — (n — 
1) coth [(3(pL — U2)\ as a function of n for various tem- 
peratures, in the non-interacting and atomic limits. It 
is clear that C is not a function of temperature only, as 
stated in Ref. ^ but varies with n and U. 

In conclusion, we have shown that the symmetry prop- 
erties (0) and (||), obtained in Refs. 0-|[ can be derived 
by means of the equation of motion and are valid for any 
temperature, including also zero temperature. These re- 
lations are exact relations and are valid for any dimen- 
sion of the system; for any value of the Hubbard inter- 
action U and of the applied magnetic field h. Further- 
more, they relate bosonic and fermionic propagators. In- 
deed, any approximation method should satisfy them in 
order to treat on equal footing one- and two- particle 
Green's functions, and to preserve spin and pseudo-spin 
symmetries^. We have also shown that the solutions pro- 
posed in Refs. [I] and || for the magnetization and for the 
particle density are not valid. 
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